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Truncation of stellar disks in galaxies at z « 1 * 
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Abstract. We report here the first evidence for stellar disk truncation at high redshift, based on surface photometry of a sample 
of 16 high redshift (0.6 < z < 1.0) disk galaxies from the GOODS HST/ACS data. The radial profiles are best fit by a double 
exponential profile. This result agrees with the profile of disks in local galaxies. The cosmological surface brightness dimming 
at this redshift range only allows us to detect galaxies with spatially "early" truncation, Rj, r /hj n ^ 3.5. Six galaxies show the 
radial double exponential structure, with an average value of Rb r /hi„ ~ 1.8. Such "early" truncated galaxies are missing in local 
samples so far. This result opens the ground for observing directly disk evolution through the study of the truncation radius as 
a function of redshift. 
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1. Introduction 

Th e finite siz e of stellar disks was first reported 
by Ivan der Kruitl dl979l) who showed that disks were sharply 
truncated at around 4 scale-lengths. Recently, evidence has 
arisen that this truncation does not occur sharply but instead 
radial profiles are best fit by a doubl e exponential shape 
(de Griis et al ., 2001; Pohle n et all 120021) . For nearby objects, 
the profiles obtained from a sample of edge-on galaxies 
show a two-slope expone ntial with Rbr/h;,, in the range 
3-5 dKregel et all 120021) . where Rt, r , the break radius, is 
the radius where the slope changes and h; n is the inner disk 
scale-length (see Pohlen et al., 2004; for a recent review on the 
truncation radius). This truncation might b e related to a surface 
density threshold for star formation jKennicutdH989l) or might 
be related to the maximum specific angular momentu m of th e 
sphere from which the disk collapsed (Ivan der Kruitl 11987). 
While galaxy formation models have b een able to reproduce 
the observed double e xponential profile (Ivan den BoschlEoOll: 
IZhang & WvseL EoO O). none of these models have considered 
evolution of the truncation radius with redshift. At z * 1, 
galaxies are half the age of local objects and one can probe 
into a period wh ere evolutionary effects may be significant 
(Ma o et all fl998 ). A study of the evolution of the truncation 
radius can help to discern between truncation arising as a relic 
from galaxy formation or due to a star formation threshold. 
Such a project has been carried out as a pilot survey using HST 
archival data and is reported in this Letter. 
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* Observations carried out using the Very Large Telescope at the 
ESO Paranal Observatory, Program ID 170.A-0788 



2. Sample and profile fitting 

The cosmological dimming (oc(l+z) 4 ) determined our choice 
of the maximum redshift to reach the outermost isophotes 
(fis ~ 25.5mag. arcsec~ 2 ). The lower redshift limit was cho- 
sen to ensure a reasonable number of objects in the final 
sample, at a reds hift where evo lutionary effects should still 
be significant iMao et all ([l998). A range of 0.6 < z < 1.0 
was finally chosen. This range implies a dimming between 
2-3 magnitudes with an isophotal threshold range between 
25.5 - 24.5 mag. arcsec~ 2 . Our sample consists of 16 galax- 
ies selected from the CDF-S and the HDF-N from the VI . 
GOODS HST/ACS Data Release JGiavalisco & et all |2004 . 
The HST/ACS z-band images were used for the analysis to en- 
sure that the rest frame of the above galaxies (V and fi-band) 
corresponds to the wavelength of widely observed local sam- 
ples. The redshifts for the southe rn fields we re obtained from 
the ESO/V LT GOODS JVanzella et alll2004T) and the COMBO 
17 projects (Wol f et alll20 04) and for the northern fields, from 
Keck DEIMOS spectra dWirth et alll2004 . The galaxies were 
visually selected to be relatively symmetric and not strongly in- 
teracting. Hubble types range from Sa to Sc (our visual classi- 
fication of the HST/ACS z-band images) and galaxy inclination 
angles (derived from the ellipticity of the outermost isophote) 
are all less than 35°. It should be emphasised that the sam- 
ple was not created with the intention of being used as a com- 
plete statistical sample. It is, in fact, an initial sample to deter- 
mine whether we can find stellar disk truncation at high-z and 
whether the derived values fall within the local range. 

Ellipse fitting to the light distribution was performed using 
the ELLIPSE task within IRAF. All the galaxies showed point- 
like nuclear regions. So, the center was fixed using the coordi- 
nates obtained by fitting a Gaussian to the nucleus. The position 
angle (PA) and the ellipticity (<?) were left as free parameters in 
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the fitting. Profiles with fixed PA and e were also obtained, no 
influence in the break radius was found. The free PA-e fitting 
profiles were chosen to determine whether there were any sys- 
tematic relation between the position of morphological features 
with the truncation radius. No attempt was made to do a bulge- 
disk decomposition, and the disk was assumed to dominate the 
light profile beyon d » 1 scale-len gth, a good assumption for 
Sb and later types dde JongMl996l) : implications of this fitting 
procedure are discussed in Section[3] The zero point to convert 
the count rate into an AB magnitude for the GOODS z-band 
was taken from the ACS web page (F850LP = 24.84). The sky 
background was subtracted by taking the mean value of boxes 
placed around the galaxies in regions far from the target galax- 
ies. The standard deviation of the median of the distribution in 
the different boxes was adopted as the error in the background 
determination. 

We then fit two exponential functions of the form 
/(R) in/out _ /^/° ut eX p(-R/h in/out ) to the 1-Dazimuthally av- 
eraged profile (in refers to the inner disk and out to the outer 
exponential disk). The criterion used for detection of trunca- 
tion was that the outer exponential should extend to at least 2 
x h OH , before reaching the noise level. The break radius was 
first visually determined, with an error around 12%. This value 
was used to separate the fitting regions for the inner and the 
outer exponentials; excluding the bulge region. The final break 
radius is defined as the radius where the two fitted exponential 
profiles cross. Fig.l shows the fitted functions plotted on top of 
the derived profiles. While the error in the inner scale-length 
is dominated by the effect of the bulge, the outer scale-length 
error is dominated by the determination of the break radius and 
the sky-subtraction. Table 1 shows the errors in the outer scale- 
length determination when doing the fitting to profiles ± 1 <x of 
the sky-level including the error in the break radius determina- 
tion. 

A correction for cosmological dimming was applied 
(H = 71 km s~ l Mpc-\ Q M = 0.3 and Q A = 0.7 through- 
out this work including converting the measured angular sizes 
to physical sizes using the gi ven redshifts), with k-c orrections 
estimated using the results of Fukugit a et alJ (11995). First, the 
k-correction for the HST/ACS z-band was computed for galaxy 
morphological type Sbc for each of the redshift values, the error 
in the K-correction is ±0.25 mag. Then, the z-band was trans- 
formed to 2?-ban d assuming the rest-frame (B - z) = 2.14+0.15 
jFukudtaet allll995l) for the same morphological type. 

3. Results and comparison to nearby objects 

Of the 16 galaxies, six show a two-slope profile that we iden- 
tify with a truncated disk. This result agrees with that found for 
nearby objects JPohlenlEoOll and references therein). The list 
of galaxies and the 1-D exponential fitting parameters and er- 
rors are presented in Table 1 . The azimuthally averaged (360°) 
profiles are shown in Fig. 1. Averaged profiles of different 
image segments were also computed (90° and 180° sections) 
to check for systematic differences but no significant changes 
were found in the profiles. The Rbr/hin occurs at 1.8± 0.5. 
Unsharp-masking was used to locate the position of the end 
of the spiral structure and to compare it to the Rb, radius. In 



two cases (123637+621 159 and 033251-27504) the truncation 
seems to be related to the end of the spiral arms. The surface 
brightness at which the break occurs happens at a mean value 
of 21.7 mag. arcsec~ 2 (rest-frame B-band). 

To compare the value obtained for the high-z galaxies 
(Rbi/hi n ranging between 1.3-2.2) with the value for nearby ob- 
jects we used two different local samples. The mean value of 
Rbi/hin obtained for a s ample of three nearby face-on galax- 
ies JPohlen et all 120021) is 3.9± 0.7. Another result from the 
analysis is that the hj n /h out ratio obtained for the high-z galax- 
ies (2_^_0;9Jj^rejriarkably similar (2 + 0.2) to that obtained 
bv IPohlen et alJ J2002I) . Notice the small average Rbr/hin ob- 
tained for the high-z galaxies compared to the value for local 
galaxies; this will be referred as "early" truncation throughout 
the paper. It is interesting to compare to this small sample be- 
cause of the similarities to our sample in the methodology used 
and the galaxy characteristics. A larger sa mple of nearby edge- 
on galaxies (20 galaxies) was studied by Kreg el et al 
with a mean value of Rbi-/hi n = 4.0± 1 . 1 . To make a fair compar- 
ison we have t o convert thes e values, obtained in the /-band, to 
the B-band lldeC.riisl [l998). The ratio h; n B/hi n; i has been esti- 
mated to be 1.32 ± 0.24 for all morphological types dde GriisL 
Il998l) . taking into account this value we can re-calculate the lo- 
cal value for the fi-band, giving a average Rbi/h; n = 2.9+0.9. In 
this way the highest values found for the high-z galaxies would 
fall in the low-value tail of the local galaxies; however, the low- 
est values still fall below the local distributi on. One should a lso 
point out that the methodology used by K regel et alJ ( 120021) to 
extract the different parameters is different to the one adopted 
in this work. Another point to notice is that there is still no clear 
result about the dependence of the break radius with wave- 
length in local galaxies. 

The remaining ten galaxies show no truncation down to 
Rb r = 3.5 - 4 h (depending on the galaxy). The derived 
scale-lengths and central surface brightnesses from the one- 
exponential fitting are comparable to those of the "early" trun- 
cated galaxies. These galaxies either do not show truncation or 
might be the counterparts of local galaxies, with spatially "late" 
truncation. 

The central disk surface brightness derived for all 16 galax- 
ies of the sample is a few magnitudes bright er than the Free man 
value 21.6 -Ji B res t-frame - 2-0 +1.0 mag ( Freemarl ll970i). No 
bulge-disk decomposition was carried out, this can lead us to 
overestimate the obtained central surface brightness. From a 
test performed on a sample of local galaxies ranging from Sa 
to Sc we estimated the maximum error from not taking into 
account the bulge to be « 0.2 mag. The errors in our central 
surface brightness estimates are high due to the k-correction 
and the bulge-disk decomposition errors. However, correcting 
for this still cannot account for the difference of around 2 mag- 
nitudes found for the central surface brightness. 

As an internal check, we determ ined the same parameters 
for the three face-on galaxies of IPohlen et alJ(Eo0 2) using their 
images and a comparable strategy to the one used for the high- 
z galaxies. We obtained similar parameters for the break ra- 
dius a nd the scale-lengths as the ones derived by IPohlen et alJ 
(l2002h . although there are some differences in the techniques 
used. 



Table 1. Galaxy parameters 
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ID 


z 


h in 




error 


hout 


error 




error 


Vo-B rl „r„, ml . 


Mb 

^» restftame 


Rbr/hin 


error 






(kpc) 


(arcsec) 


(arcsec) 


(arsec) 


(arcsec) 


(arcsec) 


(arcsec) 








123610+621334 


0.69 


6.1 


0.59 


±0.14 


0.33 


±0.10 


0.99 


±0.18 


19.7 


21.7 


1.7 


±0.5 


123637+621159 


0.78 


11.0 


0.99 


±0.13 


0.36 


±0.15 


1.32 


±0.08 


20.6 


21.9 


1.3 


±0.2 


033251-275044 


0.98 


3.5 


0.28 


±0.04 


0.15 


±0.07 


0.51 


±0.07 


18.9 


20.7 


1.8 


±0.3 


033233-274410 


0.67 


4.6 


0.45 


±0.2 


0.24 


±0.09 


0.97 


±0.05 


19.6 


21.8 


2.2 


± 1.0 


123709+622006 


1.01 


4.5 


0.36 


±0.06 


0.18 


±0.08 


0.48 


±0.08 


19.2 


20.5 


1.3 


±0.3 


123708+621252 


0.84 


5.9 


0.51 


±0.03 


0.25 


±0.13 


1.10 


±0.16 


19.6 


21.7 


2.2 


±0.3 



4. Discussion and summary 

High-z galaxies are dustier than local galaxies. This will in- 
fluence the hi n making it appear larger. The true hj n would 
need to be 50% smaller in order to raise the minimum Rbr/hin 
for the high- z galaxies to the lowest values seen locally, 
Rbr/hi n -2.6 dKregel et all [2002) (« 2 in the fi-band). In the 
nearby universe, the scale-lengths sho w a systematic increas e 
at shorter wave-lengths but not steeply (Bec kman et allfT996l) . 
It is not clear what is the dust contribution and distribution 
in galaxies at high-z which prevent us at this point to predict 
which would be the expected Rb r /hi n for the high-z galaxies. 
It is interesting to notice that the ratio between the inner and 
the outer scale-length seems to remain unchanged for the high- 
z sample which suggests that the dust does not make Rbr/hin 
appear smaller at higher redshifts. 

In a recent paper, Pohl en et alJ (2004) report a connection 
between "early" truncation and bars, where the truncation is 
associated with the outer Lindblad resonance (OLR). This im- 
plies that the break happens at Rbi/hj n » 1 — 2, consistent 
with the values found here. However, there is no clear evidence 
for the presence of bars in the 6 galaxies of our sample. The 
isophotes for the bar region are characterised by a constant PA 
whil e the e plot increases, r eaching a maximum at the end of the 
bar dWozniak et allll995l) . We do not observe this behaviour 
for the sample galaxies. However, dust might be obscuring the 
bar in the rest B-band. We could then also estimate, using the 
relation between the OLR and the truncation and relating the 
position of the OLR to the galaxy morphology, whether the 
systems are fast or slow rotators. One should also point out 
that some clearly barred objects did not show s ign of trunca- 
tion. T his fact is also observed in nearby objects (tPohlen et all 
2004) and should be further investigated. 

In summary, a double exponential provides a good fit to the 
profiles, similar to nearby galaxies. Six of the sixteen high-z 
galaxies show truncation with an average Rbi/hj n ~ 1.8, while 
for nearby objects Rbr/h; n ~ 4, reflecting several biases; our 
Rbr/hi n detection limit and the effect of dust. However, there 
are no local counterparts to the "earliest" truncated disks. The 
existence of these "early" truncated galaxies at high-z suggests 
a critical star formation density model for the origin of trun- 
cation, with disks forming inside-out. Longer wavelength ob- 
servations are needed to minimise and characterise the effects 
of dust. A larger and deeper survey could probe disk profiles 
to Rbr > 4 hin (where the break radius is observed for nearby 



objects) to obtain proper statistics on the frequency and distri- 
bution of truncation at high-z and to probe whether the 10 re- 
maining galaxies are untruncated or whether they show break 
radii similar to those of local galaxies. 
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Fig. 1. Azimuthally averaged surface brightness profiles and z-band images showing the fitted exponential functions (K-correction 
and cosmological dimming correction applied) for the 6 galaxies that show truncation. The location of the break radius has been 
overlaid on the z-band images. 



